function result=Fig1_e()
filename = 'gbsweep_199.GBs';
pathname = 'raw data/';

[data, ~] = load_GBs(filename, pathname);

data.gate = data.gate - (-20e-3); %new Vb = original Vb - Vm
n_bias = length(data.bias);
n_gate = length(data.gate);
peak = zeros(1,n_gate);
width = zeros(1,n_gate);
result.gate_2 = [];
result.peak_2 = [];
result.gate_1 = [];
result.peak_1 = [];
result.valid = [];
width_2 = [];

for i=1:n_gate
    try
        tmp = fit(data.bias',data.map(:,i),'gauss1');
    catch
        tmp.b1 = 0; tmp.c1 = 1;
    end
    width(i) = tmp.c1;
    peak(i) = tmp.b1;
    
    result.gate_2 = [result.gate_2 data.gate(i)];
    result.peak_2 = [result.peak_2 peak(i)];
    
    %remove peaks that is wider than 10mV (the peak is out of range)
    if width(i)< 10e-3 
        result.valid = [result.valid 1];
        result.gate_1 = [result.gate_1 data.gate(i)];
        result.peak_1 = [result.peak_1 peak(i)];
    else
        result.valid = [result.valid 0];
    end
    
    
end

% % 
Fcolor='k';
Bcolor='w';
Fsize=18;
Bsize=14;

f=figure();
plot(result.gate_1,-1*(1e3*(result.peak_1)+100),'LineWidth',2);

f.Color=Bcolor;
ax = gca;
ax.FontSize=Fsize;
ax.Color = 'none';
ax.XColor=Fcolor;
ax.YColor=Fcolor;
yl=ylabel('u (meV)');
xl=xlabel('Vbottom (V)');
yl.FontSize=Fsize;
xl.FontSize=Fsize;

x_limit=xlim;
x_limit=[-2.3,-0.1];
y_limit=ylim;
y_limit=[-7,20];
axis([x_limit(1),x_limit(2),y_limit(1),y_limit(2)]);


end